OTIC  FILE  COPY  ad  All  3239 


f\(LD  /bO<H.5-£L 


LIMITING  PERFORMANCE  OF  NONLINEAR 
SYSTEMS  WITH  APPLICATIONS  TO 
HELICOPTER  VIBRATION  CONTROL  PROBLEMS 


FINAL  REPORT 


WALTER  D.  PILKEY 
1982 


U.S.  Array  Research  Office 
Grant  No.  DAAG29  79  G  0015 


University  of  Virginia 


APPROVED  FOR  PUBLIC  RELEASE; 
DISTRIBUTION  UNLIMITED. 


DT!C 

|ELE'C 
APR  121982 


82  04  12  027 


UNCLASS  1 r I  ED _ 

I  SCCUWlTV  CLASSIFICATION  OF  THIS  PAGE  (*»mt  Of  gntarad) 


REPORT  DOCUMENTATION  PAGE 

READ  rNSTRUCTIONS 

BEFORE  COMPLETING  FORM 

»  REPORT NUMRCR 

UVA/ 525098  MAE  82/10 

2.  GOVT  ACCESSION  NO. 

HD-All JJIi 

4.  title  ftna  Submit) 

Limiting  Performance  of  Nonlinear  Systems  with 
Applications  to  Helicopter  Vibration  Control 
Problems 

S  type  or  REPORT  A  PERIOD  COVERED 

FINAL 

8  PERFORMING  ORO  REPORT  NUMBER 

DAG29  79  G  0015 

7.  Author^ 

Walter  D.  Pilkey 

a.  CONTRACT  or  GRANT- NUMBER!  •) 

»•  PERFORMING  ORGANIZATION  NAME  AND  AOOAESS 

<0  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  «  WORK  UNIT  NUMBERS 

ii.  Controlling  office  name  and  aooress 

U.  S.  Army  Research  Office 

Post  Office  Box  12211 

Research  Triangle  Park,  NC  27709 

12.  REPORT  DATE 

3982 

'V  NUMBER  of  PAGES 

(4.  mom i TONING  AGENCY  name  *  ADDRESS*'!/  dlltoront  from  Controlling  Ottlco) 

is.  security  class,  (oi  mu  upon) 

Unclassified 

IS*.  DECLASSIFICATION/  DOWNGRADING 
SCHEDULE 

N/A 

l«.  DISTRIBUTION  STATEMENT  (al  Jilt  Hapart) 

Approved  for  public  release;  distribution  unlimited. 

17.  distribution  statement  (al  thm  Muriel  aiiirtd  In  Slock  JO,  II  Ulfiroii  from  Rtport) 

a  a 

is.  supplementary  notes  . 

The  view,  opinions,  and/or  findir 
author(s)  and  should  not  be  cons 
position,  policy,  or  decision,  ur 

igs  contained  in 
trued  as  an  offi 
iless  so  designa 

(amrtry  oy  oioc*  nuwbmr) 

this  report  are  those  of  the 
cial  Department  of  the  Army 
ted  by  other  documentation. 

ff.  REV  WOffOS  (Continue  on  tovoemo  mi  do  H  nmcmmmory  ont 

Helicopters 

Vibrations 

Controls 

^Limiting  performance 

ABSTRACT  /CrSm  om  mooormm  •  It  nmooooomy  arf  Idontity  by  block  numbor) 

This  report  summarizes  the  accomplishments  of  a  study  exploring  new 
methods  for  the  vibration  control  of  helicopters.  Reanalysis 
methodology  permits  a  variety  of  vibration  control  problems  to  be 
solved  efficiently.  Both  analytical  and  experimental  studies  have 
been  conducted. 

1 

00,^7,1473  eornoR  or  >  mov  as  is  omoletc  UNCLASSIFIED 


^rrnmry  ri  r>r  Twit  p«r,r  rWkm »■  Pnr«*M\ 


TABLE  OF  CONTENTS 


Problem  State  and  Accomplishments  . 1 

A.  Technical  Summary  of  Analytical  Work  . 1 

B.  Technical  Summary  of  Experimental  Work  . 7 

Participating  Personnel . 21 

Publications  . 22 

Distribution  List . 23 


iDist 


PROBLEM  STATEMENT  AMD  ACCOMPLI SHMEBTS 


The  goal  of  this  program  is  to  develop  the  technology  for  the 
vibration  control  of  dynamic  systems.  Particular  emphasis  is  given 
to  helicopter  vibrations.  There  were  several  accomplishments.  A 
method  for  designing  dynamic  vibration  absorbers  to  obtain 
antiresonances  at  specified  locations  on  a  harmonically  forced 
vibratory  system  has  been  developed.  An  efficient  numerical  method 
for  the  computation  of  the  frequency  response  of  a  vibrating 
structure  as  a  function  of  its  structural  parameters  has  been  used 
to  analyze  and  design  structures  for  minimum  vibration.  A 
derivative-based  sensitivity  analysis  has  been  used  to  determine 
which  structural  element  changes  are  optimal  for  vibration 
reduction.  A  technique  for  the  optimal  reduction  of  the  steady 
state  vibration  response  over  a  frequency  range  has  been  formulated. 
The  validity  and  usefulness  of  the  methods  described  above  have  been 
justified  by  results  obtained  for  a  44  degree— of-freedom,  elastic 
line  helicopter  model.  An  experimental  helicopter  model  has  been 
built  to  investigate  active  vibration  control. 

Technical  Summary  of  Analytical  Work 

I.  Dvnamic  Vibration  Absorbers s  Antiresonance 


A  general  method  has  been  developed  for  designing  dynamic 
vibration  absorbers  which,  when  attached  to  a  harmonically  forced 
vibratory  structure,  cause  antiresonances  at  specified  points  on  the 
original  structure.  The  method  is  formulated  for  a  conservative 
multi-degree-of-freedom  system  described  by  linear  differential 
equations.  These  appendant  dynamic  systems  may  be  attached  to  the 
main  system  at  any  number  of  points.  The  formulation  is  in  terms  of 
the  receptance  of  the  absorber  at  the  attachment  point > 
consequently,  the  absorber  may  be  any  mechanical  system  whose 
receptance  is  Known  in  terms  of  the  design  parameters. 

To  describe  the  method  it  is  necessary  to  define  three  integer 
sets.  X.  A.  and  N .  X  denotes  the  set  of  degrees  of  freedom  at  which 
the  designer  wishes  to  have  antiresonance.*  A  denotes  the  set  of 
degrees  of  freedom  Where  absorbers  are  to  be  attached!  and  N  denotes 
the  set  of  points  where  external  forces  are  applied  to  the  system. 
Suppose  the  steady-state  equations  of  motion  for  the  system  in  the 
frequency  domain  are 

(K  -  w2M)  u  -  F  (11 


The  receptance  matrix  H  for  a  system  governed  by  Eq.  (l)  is  defined 
as 


H  -  (K 


-1 

(J  Ml 


(2) 
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When  the  original  system  is  considered  as  a  free  body  with 
externally  applied  forces  P^,  g?w :  and  with  the  forces  fg.  s«FA, 
exerted  by  the  absorbers  on9the  original  structure,  the  generalized 
displacement  at  each  dearee-of-freedom  iei  is  given  by 

U-EHF+EHf  f  3  > 

i  a €N  ia  g  s€A  is  s 


The  displacement  u&  at  each  degree  of  freedom  a€A  is 


u-THF  +  THf  (4  > 

a  gcN  ag  g  seA  as  s 


Let  a  denote  the  receptance  of  the  absorber  to  be  designed  and 
attached  at  degree  of  freedom  a«A.  Then 


where  f  is  the  force  exerted  by  absorber  a  on  the  original  system. 

cl 

and  is  a  function  of  the  properties  of  the  absorber  and  the 
forcing  frequency  u. 

Since  the  design  goal  is  to  maXe  each  ir  ,  i€X.  antiresonant, 
i.e..  to  have  ui  -  o  for  all  iei,  it  follows  from  Eq.  (3>  that 


EH  f  -  -  E  H  F  (6) 

seA  is  s  geN  ia  a 


Suppose  the  sets  I  and  A  both  contain  n  elements.  Equation  (6> 
represents  n  linear  equations  in  n  variables  ffl.  If  Eq.  (6}  is 

solved  for  f  ,  s?A.  and  the  solution  is  substituted  into  Eqs.  (4) 
and  ^5).  the  equation  for  is 

1 

a- - (  E  H  F  +  E  H  f  )  (7) 

a  f  aeN  ag  g  seA  as  s 

a 

Equations  (6)  and  (7)  are  the  basis  for  the  design  procedure.  For 
a  given  excitation  frequency  w.  the  forces  fg  are  found  from  Eq.  (6\ 
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and  the  receptance  is  calculated  using  Eq.  (7).  The  remaining 
problem  is  to  determine  the  structural  parameters  on  which 
depends,  such  that  a  takes  on  the  value  required  by  Eq.  (7).  These 
parameters  define  the  absorber  and  the  design  is  complete. 

II.  The  Effect  of  Structural  Properties  on  Helicopter  Response 

An  efficient  numerical  method  has  been  developed  for  the 
computation  of  the  helicopter  frequency  response  as  a  function  of 
structural  properties.  This  is  a  local  modification  method  in  which 
the  entire  helicopter  is  analyzed  once,  and  the  response  as  a 
function  of  mass  and  stiffness  parameters  is  found  by  solving  linear 
algebraic  systems  of  equations  of  reduced  order.  This  analysis 
allows  the  designer  to  reduce  vibration  levels  by  choosing  values  of 
system  property  parameters  through  a  design-by-analysis  approach 
since  the  variation  of  the  system  response  as  these  parameters  are 
changed  can  be  quickly  calculated  and  examined. 

An  outline  of  this  method  for  a  helicopter  model  made  up  of 
beam  elements  with  4x4  stiffness,  mass,  and  damping  matrices 
is  given  as  an  example.  The  receptance  matrix  R  for  a  damped 
vibratory  system  is 


R  -  (-<j2M  +  iu  C  +  K)  1  (8) 

Initially,  a  representative  set  of  system  parameter  values  is 
chosen  to  define  a  base  receptance  matrix  R.  For  the  helicopter 
described  by  this  base  receptance  matrix  R 

R-1  X  -  f  (9) 

o  o 

where  x  is  the  response  vector  and  f  the  force  vector.  Suppose 
now  tha?  the  stiffness  and  mass  matrices  of  one  fuselage  beam 
element  are  changed.  The  modified  system  is  described  by 

R*1  V  +  Dv  -  f  (10) 

•  O  "  o  o 

where  v  is  the  modified  svstem  resoonse  and  0  the  matrix  that  is 
■°  -1 

to  be  added  to  R  to  account  for  the  mass  and  stiffness  changes  in 
the  beam  element.  Equation  ( l<n  may  be  rewritten  as 

v  +  RDv  -  Rf  (11) 

•  o  ■  o  o 

The  rest  of  the  procedure  exploits  the  sparsity  of  the  matrix 
0.  Since  only  one  beam  element  has  been  changed,  the  nonzero 
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entries  of  0  form  a  4  x  4  submatrix  at  the  locations  where  the  bean 
element  contributes  to  the  global  mass  and  stiffness  matrices. 

Hence,  the  nonzero  entries  of  HD  form  an  N  x  4  submatrix  where  N  is 
the  number  of  degrees  of  freedom  of  the  system.  Consequently,  the 
product  RDy  involves  only  four  coordinates  of  the  full  response 
vector  y  .  °Let  the  4-vector  u  denote  these  four  coordinates.  Then 
it  is  possible  to  extract  from  Eq .  ( ll )  a  set  of  four  linear 
equations  to  be  solved  for  u 

u  +  Gu  -  w  ( 12 ) 

where  G  is  a  4  x  4  matrix  and  w  a  4-vector  of  known  entires.  Now, 

since  the  product  RDy  involves  only  the  entries  of  vq  given  by  u, 

it  is  clear  that,  once  u  is  solved  by  using  Eq.  (12>,  the  remaining 

elements  of  v  can  be  calculated  usinq  Eq.  (11>.  This  is  the 
-  o 

frequency  response  of  the  modified  system. 

III.  Antiresonance  using  structural  Modification 

A  structural  design  method  was  developed  to  create  antiresonances 
in  a  sinusoidally  loaded  vibrating  system.  In  this  method,  a 
structural  parameter  of  a  beam  element  in  the  helicopter  model  is 
left  as  a  design  variable  to  be  determined  such  that  at  a  specified 
degree  of  freedom  on  the  fuselage,  an  antiresonance  is  produced. 

Related  to  this  design  procedure  is  the  direct  design  of  appendant 
structures . 

Suppose  again  that  one  beam  element  is  to  be  modified.  Because 
of  the  form  of  element  property  matrices,  the  matrix  D  in  Eq.  flX) 
may  be  written  as 

D  -  «D  (13") 

o 

where  D  is  evaluated  at  some  arbitrarily  chosen  property  values, 
and  «  can  be  regarded  as  a  percent  change  in  the  structural 
properties  represented  by  the  matrix  D .  Now  Eq .  < 11 >  becomes 

V  +  <*RDv  -  Rf  ( 14 ) 

•  o  *  o  o 

and  Eq .  ( 12 't  takes  the  form 

u  +  crGu  -  w  ( 15 ) 


30  that  u  is  given  by 

u  -  ( I  +  aGI  ( 16  > 
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where  I  is  the  4x4  identity  matrix.  The  vector  u  can  be  written 
as  the  ratio  o£  two  polynomials  in  a  if  Leverrier's  Algorithm  is 
applied  to  Eq.  < 16 ) 


sh  (s) 

u .  =  — J - 

J  A(s) 


1  <  j  <  4  ( 17 ) 


where  h^(s^  is  a  third  degree  polynomial,  s  -  1/a,  and  d(s)  is  a 
fourth  degree  polynomial .  Since  the  product  aRDyQ  involves  only  the 
four  coordinates  of  y^  given  by  the  u^.  Eq.  (14)  will  yield  the 
entire  response  vector  y  as  ratios  of  polynomials 


(yo5J 


8i(s) 

A(s) 


1  <  1  <  N  (18) 


If  ( )^  is  to  be  antiresoant.  then  the  required  value  of  a  -  1/s 
is  found  by  solving  for  the  roots  of  the  polynomial  g^( s ) .  i . e . ,  for 
antiresonance 


a.  (  s )  -  0  (19) 

-k 

provided  the  root  s  is  real  and  physically  realizable.  Equation 
( 18  >  may  also  be  used  to  study  the  frequency  response  as  a  function 
of  the  design  parameter  a. 

IV.  vibration  Reduction  Over  a  Ranoe  of  Excitation  Prequencies 


As  rotor  speed  varies,  it  is  sometimes  necessary  to  consider  a 
range  of  excitation  frequencies.  Since  realistic  finite  element 
models  of  helicpters  usually  have  a  large  number  of  degees  of 
freedom,  it  becomes  important  to  be  able  to  carry  out  dynamic 
analyses  quickly  and  efficiently  as  design  variables  are  modified. 
Thus,  a  slightly  extended  version  of  the  previously  described 
reanalysis  procedure  has  been  coupled  to  a  nonlinear  programming 
optimization  method  to  arrive  at  optimal,  passive,  frequency- 
response  control  over  a  selected  interval  of  excitation  frequencies. 

A  possible  choice  for  an  objective  function  is 


0(b) 


max  ju  | 
jeJ 


(20) 


'1 
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where  the  frequency  range  of  concern  is  J  is  the  integer 

set  of  degrees  of  freedom  for  which  a  response1  reduction  is  sought ,  and 
b  is  a  vector  of  design  variables  such  as  flexural  rigidity  or 
damping  constant.  The  design  variables  are  bounded  by  physical 
realizability  criteria  so  that  the  minimization  of  the  cost  function 
G  is  subject  to  the  constraints 

<  b  <  b’J  1  <  i  <  n  (211 

along  with  the  frequency  domain  equations  of  motion 

(-w2M  +  iwC  +  R1  u  -  f  (221 


at  every  frequency  u  in  the  interval  considered. 

For  a  given  impedance  change  AZ 

AZ  -  -u2AM  +  iwAC  +  AK  {231 

Eq.  {221  may  be  rewritten  in  the  form 

[I„  +  Z”1(w1  AZ(«1]  u  -  Z"1(w1  f  (241 

where  z  is  the  imoedance  matrix  of  the  oriainal  svstem  and  I  is 
o  N 

the  Nth  order  identity  matrix.  By  talcing  advantage  of  the  sparsity 

of  the  modification  matrix  AZ.  an  exact  reanalysis  procedure 

suitable  for  automatic  machine  computation  is  derived  using  Eq. 

{241  such  that  at  any  given  iterative  step  in  the  optimization 

algorithm  is  a  reduced  order  system  of  linear  equations  which  is  solved 

instead  of  the  entire  set  (241.  The  reduced  system  has  the  form 

(I  +  D  1u,  -  u  (251 

p  p  -A  -o^ 

where  I  and  D  are  matrices  of  order  p  <  N,  the  vector  u,  is  a 
o  p  -a 

subvector  of  u.  and  the  zero  subscript  on  the  right-hand  side 

denotes  the  unmodified  system  solution.  The  elements  of  the 
response  vector  u  not  appearing  in  Eq .  (251  are  expressible  in 
terms  of  by  a  relationship  of  the  form 

u_  -  u  +  Hu.  (261 

-B  -0B  -A 

This  reanalysis  algorithm  makes  feasible  the  repeated 
evaluation  of  the  objective  function  (201  for  systems  with  large 
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numbers  of  degrees  of  freedom.  In  fact,  the  same  reanalysis 
procedure  is  also  applicable  to  the  computation  of  derivatives  with 
respect  to  design  variables  since  the  relationship 


(I 


N 


2"1(w) 

o 


aZ(w) } 


iu 

3b 


o 


3AZ(tj) 

3b 


u 


(27) 


has i the  same  form  as  Eq.  (24)  as  far  as  the  sparsity  of  the  product 
Z  (<i>)  &z(u>)  is  concerned.  The  efficient  computation  of  gradients 
of  response  with  respect  to  design  variables  is  often  of  great 
importance  in  the  use  of  format  optimization  algorithms. 

The  overall  optimization  algorithm  may  be  summarized  by  the 
following  steps: 

1)  Guess  an  initial  value  for  the  design  vector  b. 

2)  Compute  the  displacement  vector  u(x.«k)  for 

all  frequencies  of  interest  using  the  reanalysis 
procedure  given  above  in  Eqs.  (24-27). 

3)  Compute  cost  function  G(b)  using  Eq.  (20) 

4)  Compute  cost  function  derivatives  using  a  reanalvsis 
procedure  analogous  to  that  used  in  Step  2. 

5)  Apply  a  constrained  optimization  algorithm  such  as  the 
method  of  feasible  directions,  enforcing  the  constraints 
(21)  on  b.  to  compute  a  solution.  The  derivatives  of 
step  5  are  available  for  use  here. 

6)  Repeat  steps  2-5  until  a  convergence  criterion  has  been 
satisfied. 

B.  Technical  Summary  of  Experimental  Work 

A  working  model  has  been  built  to  investigate  a  concept  for 
active  vibration  control  in  helicopters.  Following  this  concept, 
the  response  of  each  blade  is  used  in  the  control  equations  to 
generate  feedback  signals  to  control  actuators  on  the  blades,  since 
"cross  signals"  are  used,  the  control  equations  can  be  written  so 
that  the  system  resembles  a  gust  alleviation  system  in  a 
conventional  airplane,  in  which  the  gust  sensor  is  mounted  on  a 
probe  ahead  of  the  wing,  otherwise,  it  can  have  a  conventional 
control  logic,  or  it  can  combine  the  conventional  system  with  the 
"gust"  system. 

A  block  diagram  giving  the  important  components  is  shown  in 
Figure  l.  while  the  present  status  of  each  component  is  discussed 
below.  This  status  is  also  summarized  in  the  figure  by  the  degree 
of  shading:  full  shading  of  a  block  indicates  completion,  half 
shading  indicates  partial  completion,  and  no  shading  indicates  that 
no  work  has  been  done.  At  the  time  of  writing  this  report,  the 
actuator  control  system  has  been  "breadboarded" .  the  final  circuit 
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board  has  been  designed  and  etched,  and  actuators  are  complete.  The 
speed  control  has  not  been  perfected  because  it  does  not  appear 
necessary. 

1.  Main  DC  Power  Supply 

A  DC  power  supply  ( Figure  2 1  is  used  to  provide  current  to 
drive  the  main  rotor.  The  current  is  controlled  by  a  bank  of 
variacs  driven  by  an  electric  motor.  AC  relays  have  been  added  so 
that  this  motor  can  be  driven  from  TTL  outputs  or  from  push  buttons 
mounted  on  the  console. 

2.  DC  Motor 

A  vertically  mounted  28  volt  DC  motor  from  an  aircraft  cabin¬ 
heating  system  is  used  to  drive  the  helicopter  rotor  directly. 

3.  Tower  (Figure  3) 

A  tower,  mounted  on  a  wooden  platform,  holds  the  DC  motor  so 
that  the  helicopter  blades  are  at  a  safe  height  above  the  ground. 

The  tower  contains  an  electronic  cabinet  and  a  shelf  for  electronic 
components . 

4 .  Flow  Environment 

The  tower  is  presently  in  a  small  room  with  thick  walls 
designed  for  rotor  experiments.  Airflow  conditions  in  the  room  are 
unfavorable,  but,  according  to  safety  standards  used  by  our  rotor 
dynamics  and  centrifugal  separation  laboratories,  we  would  otherwise 
need  at  least  1/10  inch  of  mild  steel  for  protection  should  blades 
fly  off.  Ultimately,  the  tests  will  be  run  in  moving  air.  either 
behind  a  small  wind  tunnel,  or  in  a  larger  tunnel. 

5 .  Speed  Control 

A  TTL  circuit  has  been  breadboarded  to  control  the  speed  of  the 
motor  through  the  relays  on  the  main  DC  power  supply.  It  was  found 
to  be  unstable  when  using  a  once-per-revolution  signal.  The  signal, 
from  a  proximity  probe,  has  since  been  changed  to  16  per  revolution 
by  addition  of  a  sprocket-like  wheel.  However,  the  speed  control 
has  not  been  re-activated,  largely  because  the  present  system  has 
proved  to  be  quite  steady.  The  speed  control  can  be  seen  on  the 
control  console  in  Figure  2. 

6.  Model  Helicopter  Rotor  (Figure  4) 

A  heliboy  model  helicopter  rotor  was  obtained  commercially. 

This  consists  of  the  rotor  hub.  blades,  axle,  bearings,  bearing 
blocks,  swashplate.  and  cyclic  and  collective  control  rods.  The 
assembly  was  mounted  for  direct  drive  from  the  motor,  and  provisions 
were  made  for  adjustments  to  cyclic  and  collective  pitch.  The 
so-called  "rigid  rotor"  configuration  option  was  used. 

7.  strain  Gages 

Short  aluminum  blocks  were  made  and  inserted  into  the  blade 
roots.  These  have  complete  strain  gage  bridges  capable  of  reading 
root  bending  moments,  without  being  affected  by  tension  or  torsion 
loads.  They  extend  the  blade  radius  35  mm. 
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8.  36  V  Instrumentation  Power  Supply 

A  two-amp.  36  volt,  unregulated  DC  power  supply  was  built  to 
provide  power  to  the  electronics  on  the  hub. 

9.  Slip  Rings 

Two  bronze  slip  rings  were  built  and  mounted  on  the  axle. 

These  supply  the  36  volt  power  to  the  hub.  The  axle  was  grooved  to 
admit  the  two  conductors. 

10.  Strain  Gauge  Amplifier 

A  round  circuit  board .  mounted  just  above  the  blades,  contains 
regulated  5V  and  18V  power  supplies,  the  5V  supply  being  used  only 
to  polarize  the  strain  gauge  bridges.  The  signal  from  each  bridge 
is  amplified  by  two  LM  741  operational  amplifiers.  Three  20-turn 
potentiometers  are  intailed  for  adjustments  on  each  circuit.  The 
board  shows  clearly  in  Figures  7  and  8. 

11.  DC  to  PM  Conversion 

The  amplified  DC  signals  are  converted  to  FM  by  the  VCO 
portions  of  two  LM  562B  phase- loch  loops.  Center  frequencies  are 
1.0  and  1.414  MHz. 

12.  Annular  Antennae 

The  two  FM  signals  are  transmitted  by  annular  antennae  which 
have  been  made  of  circuit  board  material  from  the  patterns  shown  in 
Figures  5  and  6.  One  of  the  antennae  is  attached  to  the  hub  below 
the  blades,  the  other  is  fixed  to  the  mounting  on  the  motor.  The 
rotating  half  has  two  180°  segments  for  the  1.0  and  1.414  MHz 
signals,  these  correspond  to  a  360*  segment  on  the  fixed  half,  which 
receives  the  combined  signal.  Corresponding  360°  segments  carry  the 
signal  ground.  The  antennae  can  be  seen  mounted  in  Figures  7  and  8. 

13.  Tims-16  Signal 

The  rotating  antennae  half  also  contains  sixteen  sectors.  It 
was  intended  that  these  would  be  used  in  conjunction  with  a  pair  of 
short  segments  in  the  fixed  half  to  generate  a  revolution  times  16 
signal  from  which  the  blade  azimuth  could  be  determined  and  supplied 
to  the  TRS— 80  computer.  However,  the  present  configuration  includes 
a  sixteen  tooth  wheel  on  the  rotor  shaft  from  which  a  fixed 
proximity  probe  picks  up  the  "times-16"  signal.  Circuits  on  the 
demodulation  board  detect  the  positive  and  negative  going  edges  of 
this  signal  and  thus  convert  it  to  a  "times-32"  signal. 

14.  Card  Rack  and  Power  Supply 

The  electronic  cabinet  on  the  tower  contains  a  card  rack  with  a 
+5V.  ±15  power  supply. 

15.  FM  to  DC  Monitor 

A  card  in  the  rack  contains  permanent  circuits  which  demodulate 
the  1.0  and  1.414  MHz  strain  gauge  signals  or  the  .861  and  1.189  MHz 
control  feedback  signals  using  LM562B  phase- lock  loops.  The 
demodulated  signals  can  be  displayed  on  a  oscilloscope,  or  they  can 
be  read  directly  on  a  voltmeter. 
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16.  Interface  Board 

A  wire— wrap  board,  hinged  to  the  top  of  an  electrical  box. 
contains  the  interface  circuits  for  the  TRS-80  Model  III .  The  box 
contains  power  supplies  for  +5V  and  ±15V. 

17.  PM  to  DC  Demodulator  and  Trigonometer  Functions 

The  PM  to  DC  demodulation  circuit  on  the  interface  board 
demodulates  the  1.0  and  1.414  MHz  strain  gauge  signals.  The  two 
outputs  are  read  into  sample  and  hold  circuits  as  x^  and  x^:  these 
are  also  multiplied  by  sines  and  cosines  of  the  azimuth  angles 
discussed  below,  using  four  DAC  7523  multiplying  D/A  converters. 

The  times-32  signal  derived  from  the  proximity  probe  generates  a 
5-bit  cosine  which  is  used  to  address  a  fusible-link  PROM  containing 
a  sine  table.  By  logical  operations  on  the  5-bit  count,  the 
following  signals  are  generated: 

x  -  sin  R.  x,  -  cos  9.  x_  •  sinf  9  +  180®  x.  »  cos(  9  +  180°  * 

where  9  is  the  blade  azimuth  angle.  The  "32-times"  signal  also  sets 
a  flip/ flop  so  that  on  completing  the  cycle,  the  computer  can  be 
programmed  to  count  for  new  data  before  starting  another  cycle. 

18.  A/D  Conversion 

A  CD AS  9S2R  A/D  analyzer  reads  both  inputs  X^  and  x^.  together 
with  their  products  with  the  sines  and  cosines  of  the  azimuth  angles 

X  .  x,  ■  x  and  X  . 

2  3  5  6 

19.  Radio  Shack  TRS-80  Model  III  Computer 

This  is  a  standard  Model  III  computer,  with  two  disc  packs,  and 
is  sold  as  a  home  computer.  The  micorprocessor  is  a  Z80.  and  there 
are  48K  of  dynamic  storage.  For  test  purposes,  it  can  be  programmed 
in  BASIC,  but  for  the  present  use.  it  must  be  programmed  in  machine 
language  so  that  reasonable  speed  can  be  achieved.  The  interface 
board  is  mapped  into  input/output,  using  instructions  shown  in  Table 
1 . 


20.  Machine  Language  Program 

The  complete  control  logic  is  incorporated  in  the  computer 
program.  Two  conflicting  approaches  to  this  have  been  considered: 

(1>  In  the  traditional  control  theory  approach,  each  control 
output  Y  ( blade  position >  would  be  generated  from  the  corresponding 
strain  gauge  reading  X. 

(  2 ^  in  the  "gust  alleviator”  approach,  each  blade  would  act  as 
a  sensor  to  control  the  other,  following  blade.  This  assumes  that 
the  disturbances  are  quasi-cyclic.  while  the  first  approach  assumes 
that  they  are  completely  random. 
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In  generating  an  initial  control  program,  an  intermediate 
approach  has  been  taken,  with  each  input  generated  from  both  strain 
gauges.  The  control  equation  is 

;  }  -  tc3<x1x2x3x4x5x6}T  ♦  D 

with  the  inouts  X. .  X, .  X, .  and  X.  used  in  such  a  wav  that 
"once-per-revolutaon"  signals,  such  as  pilots  control  inputs,  do  not 
drive  the  actuators  to  the  stops.  Values  for  the  c  and  D  matrices  can 
be  inserted  through  a  BASIC  master  program.  The  constant  terms  are 
included  because  all  input  signals  represent  positive  voltages  in  the 
range  of  0  to  5V  converted  to  8-bit  digital  numbers.  Thus,  sin  o  is 
read  as  80  hexadecimal.  At  present,  a  multiplcation  subroutine  with 
an  8-bit  output  is  used.  However,  once  some  initial  testing  has  been 
completed,  it  will  be  easy  to  modify  this  program.  One  important 
consideration  is  whether  to  continue  reading  the  four  trigonometric 
functions,  suffering  the  long  conversion  delays,  or  to  generate  these 
internally. 

TABLE  l  TRS— 80  Model  XXI  Instructions  to  Interface  Board 


rUNCTION 

ASSEMBLY  CODE 

LATCH  sin  9 

OUT  /0KA 

LATCH  cos  9 

OUT  flt.A 

SAMPLE  CH  *1 

OUT  (  2  ).A 

HOLD  CH  #1 

OUT  (3).  A 

LATCH  Sin  19  + 

180*  ) 

OUT  (4^, A 

LATCH  COSf9  + 

180°  ^ 

OUT  f5^.A 

SAMPLE  CH  *2 

OUT  (  6).A 

HOLD  Ch  #2  -  X 

A 

OUTf 7). A 

OUTPUT  CH  #1  - 

OUT  f 17  > , A 

OUTPUT  CH  #2  « 

OUT  fl8>.A 

MUX  X 

OUT  (32). A 

MUX  X*  *  cos  9 

-  X 

OUT  (33>.A 

MUX  X  *  cos  9 

-  X3 

OUT  H4).A 

MUX  X* 

OUT  (36). A 

MUX  X*  *  cost  9 

+  180°  )  -  Xg 

OUT  (37KA 

MUX  X*  *  sinf  9 

+  180°  )  -  X- 

OUT  (3B).Z 

HEAD  STATUS 

© 

IN  A.  (0) 

READ  AID 

IN  A.  / 16  ) 

CLEAR  FLIP /PLOP 

IN  A.  (32) 

GENERATE  TIMING 

SIGNAL 

IN  A.  (  48  ) 

A  -  Contests  of  computers  8-bit  accumulator 
9  -  Azimuth  angle 


21.  0/A  Converters 

Output  data  Y  .  Y+.  for  the  two  blades  are  latched  in  digital 
form  on  the  interface  board,  and  converted  to  analog  form  in  two  DAC 
7523  8-bit  D/A  converters. 

22.  DC  to  PM  Conversion 

The  analog  signals  to  the  two  blades  are  converted  to  PM 
signals  at  .841  and  1.189  MHz  using  Uf  562B  phase  lock  loops.  By 
restricting  frequency  variations  to  within  the  range  92%  to  1.09%  of 
nominal  values,  cross-talk  between  the  four  channels  should  be 
eliminated.  It  should  be  noted  that  the  VCO's  on  phase- lock  loops 
do  not  produce  pure  sinusoidal  signals,  therefore  the  receivers  can 
lock  onto  harmonics  of  the  nominal  frequencies  so  that  the  system 
is  restricted  to  a  2:1  frequency  range. 

23.  PM  to  DC  Monitor 

The  same  card  can  be  used  on  the  card  rack  to  monitor  the  1.189 
and  1.682  MHz  signals  from  the  interface  board,  as  is  used  to 
monitor  the  strain  gauge  signal. 

24.  Annular  Antennae 

The  two  PM  signals  are  to  be  transmitted  back  onto  the  hub  via 
the  same  pair  of  annular  antennae  as  are  used  to  transmit  the  strain 
gauge  signals.  The  signals  share  a  common  ground  track. 

25.  PM  to  DC  Demodulation 

A  second  circuit  board,  to  be  mounted  above  the  strain  gage 
board,  will  contain  a  demodulation  circuit,  using  LM  562B  phase-lock 
loops . 

26.  Actuator  Solenoids 

circuitry  on  the  demodulator  board  will  drive  the  actuator 
solenoids  which  have  already  been  mounted  in  place  of  the  collective 
pitch  linkages.  These  solenoids  act  one-way  only,  and  have  been 
provided  with  return  springs.  Feedback  position  sensors  will  be 
incorporated  into  a  feedback  control  system,  which  includes  four 
operational  amplifiers  and  two  transitors  to  drive  each  solenoid. 
This  system  has  already  been  evaluated,  and  the  circuit  board  is 
being  assembled. 
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PUBLICATIONS 


Several  computational  design  techniques  for  the  vibration 
control  of  helicopters  have  been  developed.  The  following 
manuscripts  describe  this  work. 

1.  "Synthesis  of  Dynamic  Vibration  Absorbers."  submitted  to 
J.  of  Engineering  for  Industry. 

2.  "Helicopter  Vibration  Reduction  by  Local  Structural 
Modification . "  Accepted  for  publication  by  the  J.  of  American 
Helicopter  Society. 

3.  "Structural  Modification  to  Achieve  Antiresonance  in 
Helicopters , "  J.  of  Aircraft. 

4.  "Eigensolution  Reanalysis  of  Rotor  Dynamic  Systems  by  the 
Generalized  Receptance  Method,"  Accepted  for  publication  in  the 
International  J.  for  Numerical  Methods  in  Engineering. 

5 .  "Optimal  Band-selective  Vibratory  Response  Reduction 
Techniques  Applied  to  Helicopters."  ready  to  submit  for  publication. 

The  contents  of  the  first  three  manuscripts  were  summarized  and 
our  more  recent  work  ( paper  5 1  was  accorded  a  more  detailed 
explanation  in  this  report. 
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